A Zoology of Bell inequalities resistant to detector inefficiency 
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We derive both numerically and analytically Bell inequalities and quantum measurements that 
present enhanced resistance to detector inefficiency. In particular we describe several Bell inequalities 
which appear to be optimal with respect to inefficient detectors for small dimensionality d — 2, 3, 4 
and 2 or more measurement settings at each side. We also generalize the family of Bell inequalities 
described in Collins et all |s| to take into account the inefficiency of detectors. In addition we 
consider the possibility for pairs of entangled particles to be produced with probability less than 
one. We show that when the pair production probability is small, one must in general use different 
Bell inequalities than when the pair production probability is high. 

PACS numbers: 03.65.Ta 



I. INTRODUCTION 

A striking feature of quantum entanglement is non- 
locality. Indeed, as first shown by Bell in 1964 jj| clas- 
sical local theories cannot reproduce all the correlations 
exhibited by entangled quantum systems. This non-local 
character of entangled states is demonstrated in EPR 
experiments through the violation of Bell inequalities. 
However due to experimental imperfections and techno- 
logical limitations, Bell tests suffer from loopholes which 
allow, in principle, the experimental data to be repro- 
duced by a local realistic description. The most famous 
of these loopholes are the locality loophole and the de- 
tection loophole. Experiments carried on photons have 
closed the locality loophole |2| and recently Rowe et al 
closed the detection loophole using trapped ions || . But 
so far, 30 years since the first experiments, both loopholes 
have not been closed in a single experiment. 

The purpose of this paper is to study how one can de- 
vise new tests of non-locality able to lower the detector 
efficiency necessary to reject any local realistic hypoth- 
esis. This could be a way towards a loophole-free test 
of Bell inequalities and is important for several reasons. 
First, as quantum entanglement is the basic ingredient 
of quantum information processing, it is highly desirable 
to possess undisputable tests of its properties such as 
non-locality. Even if one is convinced (as we almost all 
are) that nature is quantum mechanical, we can imag- 
ine practical situations where it would be necessary to 
perform loophole-free tests of Bell inequalities. For ex- 
ample, suppose you buy a quantum cryptographic device 
based on Ekert protocol. The security of your crypto- 
graphic apparatus relies on the fact that you can violate 
Bell inequalities with it. But if the detectors efficien- 
cies aren't sufficiently high, the salesman can exploit it 
and sell to you a classical device that will mimic a quan- 
tum device but which will enable him to read all your 
correspondence. Other reasons to study the resistance 
of quantum tests to detector inefficiencies are connected 



to the classification of entanglement. Indeed an impor- 
tant classification of entanglement is related to quantum 
non-locality. One proposed criterion to gauge how much 
non-locality is exhibited by the quantum correlations is 
the resistance to noise. This is what motivated the se- 
ries of works (J, |j| that led to the generalization of the 
CHSH inequality to higher dimensional systems J(i| . The 
resistance to inefficient detectors is a second and different 
criterion that we analyse in this paper. It is closely re- 
lated to the amount of classical communication required 
to simulate the quantum correlations 0. 

The idea behind the detection loophole is that in the 
presence of unperfect detectors, local hidden variables 
can "mask" results in contradiction with quantum me- 
chanics by telling the detectors not to fire. This is at 
the origin of several local hidden variable models able 
to reproduce particular quantum correlations if the de- 
tector efficiencies are below some threshold value 77* (see 
H ||, [l(| for example). In this paper, we introduce two 
parameters that determine whether a detector will fire or 
not: 77, the efficiency of the detector and A, the probabil- 
ity that the pair of particles is produced by the source of 
entangled systems. This last parameter may be impor- 
tant for instance for sources involving parametric down 
conversion where A is typically less than 10%. So far, 
discussions on the detection loophole where concentrat- 
ing on 77, overlooking A. However we will show below 
that both quantities play a role in the detection loophole 
and clarify the relation between these two parameters. In 
particular we will introduce two different detector thresh- 
olds: 77*, the value above which quantum correlations ex- 
hibit non- locality for given A, and 77* , the value above 
which quantum correlations exibit non-locality for any A. 

We have written a numerical algorithm to determine 
these two thresholds for given quantum state and 
quantum measurements. We then searched for optimal 
measurements such that 77^ =1 and 77* acquire the lowest 
possible value. In the case of bipartite two dimensional 
systems the most important test of non-locality is the 
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TABLE I: Optimal threshold detector efficiency for varying 
dimension d and number of settings (N a x Nt) for the detec- 
tors. r)f =1 is the threshold efficiency for a source such that the 
pair production probability A = 1 while r^ A is the threshold 
efficiency independent of A. The column p gives the amount 
of white noise p that can be added to the entangled state so 
that it still violates locality (we use for p the same defini- 
tion as that given in P). The last column refers to the 
Bell inequality that reproduce the detection threshold. Com- 
pletely new inequalities introduced in this paper are indicated 
by "New". 



d 


N a x N b 


i^ 1 




P 


Bell inequality 


2 


2x2 


0.8284 


0.8284 


0.2929 


CHSH 


2 


3x3 


0.8165 




0.2000 


New 

(see also ref jl[ H ) 


2 


3x3 




0.8217 


0.2859 


New 


2 


3x4 




0.8216 


0.2862 


New 


2 


4x4 




0.8214 


0.2863 


New 


3 


2x2 


0.8209 


0.8209 


0.3038 


based on ref.Q 


3 


2x3 


0.8182 


0.8182 


0.2500 


New 
(related to ref 


3 


3x3 


0.8079 




0.2101 


New 


3 


3x3 




0.8146 


0.2971 


New 


4 


2x2 


0.8170 


0.8170 


0.3095 


based on ref. pi 



2x3 
3x3 



0.8093 0.2756 
0.7939 0.2625 



New 
New 



5 2 x 2 0.8146 0.8146 0.3128 based on ref.[ 

6 2 x 2 0.8130 0.8130 0.3151 based on ref.[ 

7 2 x 2 0.8119 0.8119 0.3167 based on ref.[ 
oo 2 x 2 0.8049 0.8049 0.3266 based on ref. [ 



CHSH inequality jlj}] . Quantum mechanics violates it if 
the detector efficiency r\ is above = 2/(v / 2 + 1) w 0.8284 
for the maximally entangled state of two qubits. In the 
limit of large dimensional systems and large number of 
settings, it is shown in H that the efficiency threshold 
can be arbitrarily lowered. This suggests that the way 
to devise optimal tests with respect to the resistance to 
detector inefficiencies is to increase the dimension of the 
quantum systems and the number of differents measure- 
ments performed by each party on these systems. (This 
argument will be presented in more details in Q). We 
have thus performed numerical searches for increasing 
dimensions and number of settings starting from the two 
qubit, two settings situation of the CHSH inequality. 
Our results concern "multiport beam splitters measure- 
ments" |l2] ] performed on maximally entangled states. 
They are summarized in Table I. Part of these results 
are accounted for by existing Bell inequalities, the other 
part led us to introduce new Bell inequalities. 

The main conclusions that can be drawn from this 
work are: 

1. Even in dimension 2, one can improve the resistance 



to inefficient detectors by increasing the number of 
settings. 

2. One can further increase the resistance to detection 
inefficiencies by increasing the dimension. 

3. There are different optimal measurements settings 
and Bell inequalities for a source that produces en- 
tangled particles with high probability (A k, 1) and 
one that produces them extremely rarely (A — *■ 0). 
Bell inequalities associated with this last situation 
provide a detection threshold that doesn't depend 
on the value of the pair production probability. 

4. For the measurement scenarios numerically accessi- 
ble, only small improvements in threshold detector 
efficiency are achieved. For instance the maximum 
change in threshold detector efficiency we found is 
approximatively 4% 

The paper is organized as follows: First, we review 



briefly the principle of an EPR experiment in section [I A 
and under which condition such an exp eriment admits a 



we 
to 



local-realistic description in section II B. In section HC 
we clarify the role played by ry and A in the detection 
loophole. We then present the techn ique we used to per- 
form the numerical searches in II D and to con struc t the 
new Bell inequalities presented in this paper in |IIE|. S ec- 
tion [II contains our results. In particular in |III A 
generalize the family of inequalities introduced in jp 
take into account detection inefficiencies and in [ill C we 
present the two different Bell inequalities associated to 
the two-dimensional three by three settings measurement 
scenario. In the Appendix, we collect all the measure- 
ment settings and Bell inequalities we have obtained. 



II. GENERAL FORMALISM 

A. Quantum correlations 

Let us review the principle of an a EPR experiment: 
two parties, Alice and Bob, share an entangled state 
Pab ■ We take each particle to belong to a d dimensional 
Hilbert space. The parties carry out measurements on 
their particles. Alice can choose between N a different 
von Neumann measurements A4 (i = 1, . . . , N a ) and Bob 
can choose between Nb von Neumann measurements Bj 
(j = 1, . . . , Nb). Let k and I be Alice's and Bob's out- 
comes. We suppose that the number of possible outcomes 
is the same for each party and that the values of k and 
I belong to {0, . . . , d — 1}. To each measurement Aj is 
thus associated a complete set of d orthogonal projectors 
A\ = \A^)(A^\ and similarly for B } ■. Quantum mechanics 
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predicts the following probabilities for the outcomes 



with the conditions: 



P? M {B ) = Tr((l A (8 B l j)p ab ) , 
P^iAi) = Tr((4 fc ®l B K b ). 



(1) 



In a real experiment, it can happen that the measure- 
ment gives no outcome, due to detector inefficiencies, 
losses or because the pair of entangled states has not 
been produced. To take into account these cases in the 
most general way, we enlarge the space of possible out- 
comes and add a new outcome, the "no-result outcome" , 
which we label 0. Quantum mechanics now predicts a 
modified set of correlations: 





= k, B 3 


-0 


= tfpW^Bj) k,l^ 




= Mj 


= 


= Xnii-^p^iBj) i + 




= k,B j 


= 0) 


= \ v (l- v )P? M (A t ) k + 




= <D,B j 


-0) 


= l-A + A(l-77) 2 . 



(2) 

where rj is the detector efficiency, and A is the probability 
that a pair of particles is produced by the source of en- 
tangled systems. By detection efficiency rj we mean the 
probability that the detector gives a result if a particle 
was produced, i.e. rj includes not only the "true" effi- 
ciency of the detector but also all possible losses of the 
particle on the path from the source to the detectors. 



B. Local Hidden Variable Theories &: Bell 
Inequalities 



\) of an EPR ex- 
hidden variable 



Let us now define when the results ( 
periment can be explained by a loca 
(lhv) theory. In a lhv theory, the outcome of Alice's 
measurement is determined by the setting Ai of Alice's 
measurement apparatus and by a random variable shared 
by both particles. This result should not depend on the 
setting of Bob's measurement apparatus if the measure- 
ments are carried out at spatially separated locations. 
The situation is similar for Bob's outcome. We can de- 
scribe without loss of generality such a local variable the- 
ory by a set of ( d+ 1 ) Na +Nh probabilities p Kl .. . K Na L 1 ...L Nb 
where Alice's local variables Ki s {0, ... , d — 1, 0} spec- 
ifies the result of measurement Ai and Bob's variables 
Lj e {0, ... , d — 1, 0} specify the result of measurement 
Bj. The correlations P(Ai — K,Bf. — L) are obtained 
from these joint probabilities as marginals. The quan- 
tum predictions can then be reproduced by a lhv theory 



if and only if the following N a Nb(d 
obeyed: 



1) equations are 



EWi^/^L = = K, Bj = L) (3) 

KL 



E PKL = f 
KL 

PKL > , 



(4) 
(5) 



where we have introduced the notation K = K\ . . . K^ a 
and h = Li . . . Lff b . Note that the equations (Q) are not 
all independent since quantum and classical probabilities 
share additional constraints such as the normalization 
conditions: 



P(Ai = K, Bj = L) = 1 



(6) 



K.L 



or the no-signalling conditions: 

P{A i =K)=Y,P{A i =K,B j = L) Vj (7) 



and similarly for Bj. 

An essential result is that the necessary and sufficient 
conditions for a given probability distribution pQ M to be 
reproducible by a lhv theory can be expressed, alterna- 
tively to the equations (^|) , as a set of linear inequalities 
for p ( 3 M j the Bell inequalities. They can be written as 



where 



/ — I rr + I$ r + I r + /00 < c 



Z - = E E cf j P(A i =k,B j =l) 

i,3 k,l^0 

hr = £5>$P(^ = 0,B; = O 

In = $>ffP(A i = 0,1^ = 0). 



(8) 



(9) 



For certain values of rj and A, quantum mechanics can 
violate one of the Bell inequalities (|J) of the set. Such a 
violation is the signal for experimental demonstration of 
quantum non-locality. 



C. Detector efficiency &: pair production 
probability 

For a given quantum mechanical probability distri- 
bution pQ M and given pair production probability A, 
the maximum value of the detector efficiency r\ for 
which there exists a lhv variable model will be denoted 
V*(P QM )- K has been argued §, [f| that r?* should not 
depend on A. The idea behind this argument is that 
the outcomes (0, 0) obtained when the pair of particles is 
not created are trivial and hence it seems safe to discard 
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them. A more practical reason, is that the pair produc- 
tion rate is rarely measurable in experiments. Whatever, 
the logical possibility exists that the lhv theory can ex- 
ploit the pair production rate. Indeed, we will show be- 
low that this is the case when the number of settings of 
the measurement apparatus is larger than 2. This moti- 
vates our definition of threshold detection efficiency valid 
for all values of A 



vt X = maxfo*) = lim 77* 



(10) 



The second equality follows from the fact that if a lhv 
model exists for a given value of A it also exists for a 
lower value of A. 

Let us study now the structure of the Bell expression 
I(QM) given by quantum mechanics. This will allow us 
to derive an expression for rj^ x . Inserting the quantum 
probabilities (^) into the Bell expression of eq. (||) we 
obtain 

I(QM) = \ v 2 I rr (QM) + Xr){l - n)Ifr(QM) 
+ Ar?(l - v)IMQM) + (1 + Hv 2 - 27?)) J2 4f (11) 

id 

where I®, M is obtained by replacing P(Ai — k, Bj = I) 



with (Ai,Bj) in I rr and J^, by replacing P(Ai = 
<D,B j = l) with P^ M (Bj) in 7 0r and similarly for I^ M . 

For 77 = 0, we know there exists a trivial lhv model and 
so the Bell inequalities cannot be violated. Replacing 77 
by in (fLI|) we therefore deduce that 



< c. 



This divides the set of Bell inequalities into two groups: 

those such that Y\- • cff < c and those for which 
* — iL ij l j 

J2i j c ij = c - Let us consider the first group. For small 
A, these inequalities will cease to be violated. Indeed, 
take rj = 1 (which is the maximum possible value of the 
detector efficiency), then (O) reads 



I(QM) = XlQ r M + (i-A)£\ 



(13) 



The condition for violation of the Bell inequality is 
I(QM) > c. But since c ij < c ; l° r sufficiently small 
A we will have I(QM) < c and the inequality will not be 
violated. These inequalities can therefore not be used to 
derive threshold rj^ x that do not depend on A, but they 
are still interesting and will provide a threshold 77^ de- 
pending on A. Let us now consider the inequalities such 
that J2i j c ij = c - Then A cancels in ( |TT| ) and the condi- 
tion for violation of the Bell inequality is that 77 must be 
greater than 



2c - I, 



QM 
0r 



tQM 



t QM 



',•0 



It is interesting to note that if quantum mechanics vi- 
olates a Bell inequality for perfect sources A = I and per- 
fect detectors 77 = I, then there exists a Bell inequality 
that will be violated for 77 < 1 and A — > 0. That is there 
necessarily exists a Bell inequality that is insensitive to 
the pair production probability. Indeed the violation of 
a Bell inequality in the case A = I, 7/ = I implies that 
there exists a Bell expression I rr such that I rr {QM) > c 
with c the maximum value of I rr allowed by lhv theories. 
Then let us build the following inequality 



I = I r 



Ir$ + hr + £ cffPiA, = 0, Bj = 0) < 



(15) 



where ^ c®j = c and we take in I T $ and I§ r sufficiently 



negative terms to insure that I < c. For this inequality, 



(2c — I, 



QM _ tQM n i( jQm _ tQM _ T QM x -. 
0r r0, )/\ c + 1 rr Hr V0 ' . ' 

which shows that Bell inequalities valid VA always exist. 
One can, in principle, optimize this inequality by ta king 
J r and I% T as large as possible while ensuring that ( |f5| ) 
is obeyed. 

From the experimentalist's point of view, Bell tests in- 
volving inequalities that depend on A need all events to 
be taken into account, including (0, 0) outcomes, while 
in tests involving inequalities insensitive to the pair pro- 
duction probability, it is sufficient to take into account 
events where at least one of the parties produces a re- 
sult, i.e. double non-detection events (0, 0) can be dis- 
carded. Indeed, first note that one can always use the 
normalization conditions (^|) to rewrite a Bell inequality 
such as (|^) in a form where the term Jgg does not ap- 
(12) pear. Second, if the events (0,0) are not recorded in an 
experiment, the measured probabilities are relative fren- 
quencies computed on the set of all events involving at 
least one result on one side. The probabilities measured 
in such experiments can be obtained from the probabili- 
ties (||) by replacing A with A' = A/(I-(I-77) 2 ). While 
this rescaling of A is legitimate for inequalities that do 
not depend on the value of A, it is however incorrect to 
pcrfom it for inequalities depending of A, in particular 
this will affect the detection threshold. 



D. Numerical search 

We have carried numerical searches to find measure- 
ments such that the thresholds ?7^ =1 and r% x acquire the 
lowest possible value. This search is carried out in two 
steps. First of all, for given quantum mechanical proba- 
bilities, we have determined the maximum value of 77 for 
which there exists a local hidden variable theory. Second 
we have searched over the possible measurements to find 
the minimum values of 77*. 

In order to carry out the first step, we have used the 
fact that the question of whether there are classical joint 
(14) probabilities that satisfy ^ with the conditions (|^J|) 
is a typical linear optimization problem for which there 
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exist efficient algorithms [|14| . We have written a program 
which, given A, 77 and a set of quantum measurements, 
determines whether (||) admits a solution or not. rj X is 
then determined by performing a dichotomic search on 
the maximal value of 77 so that the set of constraints is 
satisfied. 

However when searching for r^ x it is possible to dis- 
pense with the dichotomic search by using the follow- 
ing trick. First of all because all the equations in eq. 
(0) are not independent, we can remove the constraints 
which involve on the right hand side the probabilities 
P{Ai = 0,Bj = 0). Second we define rescaled variables 
A(l — (1 — ?7)^)pkl = Pkl- Inserting the quantum prob- 
abilities eq. (0) we obtain the set of equations 



^2pKLS Ki ,kSL it i = aP^ M (Ai, Bj) k, I £ 

KL 

/JPKL^Kf,0^7,' 



KL 



$>KL<^,^,0 = (1 - -)P« M (A i ) k + 



KL 



with the normalization 



KL 



A 1 - (1 -rj) 2 



(16) 



(17) 



where a = tj 2 /(l — (1 — rf) 2 ). Note that A only appears in 
the last equation. We want to find the maximum a such 
that these equations are obeyed for all A. Since < A < 1 
[fill, we can replace the last equation by the condition 



X^PKL > 1. 



(18) 



KL 



We thus are led to search for the maximum a such that 
eqs. ([l6|) are satisfied and that the pkl are positive and 
obey condition ([l8|). In this form the search for rj^ x has 
become a linear optimization problem and can be effi- 
ciently solved numerically. 

Given the two algorithms that compute rj x=1 and rj^ x 
for given settings, the last part of the program is to 
find the optimal measurements. In our search over the 
space of quantum strategies we first considered the maxi- 
mally entangled state — X^ m =o l m )al TO )& m dimension 
d. The possible measurements A4 and Bj we consid- 
ered are the "multiport beam splitters" measurements 
described in ]l2| and which have in previous numerical 
searches yielded highly non local quantum correlations 
, P . These measurements are parametrized by d phases 



b d A .) and {(f) B . , . . . <j> B . ) and involve the follow- 



ing steps: first each party acts with the phase <f>Ai (jn) 
or (f>B-{fn) on the state |m), they then both carry out a 



discrete Fourier transform. This brings the state \& to: 
^ d-l 

* = im E ex P 



k,l,m—0 



27T n n 
+—m(k — L) 



i y4> Al (m) - 4>Bj (m) 
\k) a \l) b (19) 



Alice then measures \k) a and Bob |Z);,. 
probabilities dil) thus take the form 



d-l 



P? M (A t ) = 
P? M (B 3 ) = 



-3 1 Yj e*p[i(Wr 

m=0 

+ ^( fc -0)]l 2 

l/d 
l/d 



The quantum 



0s 3 (m) 



(20) 



The search for minimal ?y* =1 and ij^ x then reduces to 
a non-linear optimization problem over Alice's and Bob's 
phases. For this, we used the "amoeba" search procedure 
with its starting point fixed by the result of a randomized 
search algorithm. 

Note that these searches are time-consuming. Indeed, 
the first part of the computation, the solution to the lin- 
ear problem, involves the optimization of (d + 1)^+^ 
parameters, the classical probabilities pkl (the situation 
is even worse for rj X , since the linear problem has to be 
solved several times while performing a dichotomic search 
for r] x ). Then when searching for the optimal measure- 
ments, the first part of the algorithm has to be performed 
for each phase settings. This results in a rapid exponen- 
tial growth of the time needed to solve the entire problem 
with the dimension and the number of settings involved. 
A second factor that complicates the search for optimal 
measurements is that, due to the relatively large number 
of parameters that the algorithm has to optimize, it can 
fail to find the global minimum and converge to a local 
minimum. This is one of the reasons why, as a first step, 
we restricted our searches to "multiport beam splitter" 
measurements since the number of parameters needed to 
describe them is much lesser than that for general Von 
Neumann measurements. 

Our results for setups our computers could handle in 
reasonable time are summarized in table I. In dimension 
2, we also performed more general searches using von 
Neuman measurements but the results we obtained where 
the same as for the multiport beam splitters described 
above 



E. Optimal Bell inequalities 

Upon finding the optimal quantum measurements and 
the corresponding values of ?y* , we have tried to find the 
Bell inequalities which yield these threshold detector ef- 
ficiencies. This is essential to confirm analytically these 
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numerical results but also in order for them to have prac- 
tical significance, ie. to be possible to implement them 
in an experiment. 

To find these inequalities, we have used the approach 
developped in ||. The first idea of this approach is to 
make use of the symmetries of the quantum probabili- 
ties and to search for Bell inequalitites which have the 
same symmetry. Thus for instance if P(Ai — k, Bj = 
I) = P(Ai = k + m mod d, Bj = I + m mod d) for all 
m G {0, . . . , d — 1}, then it is useful to introduce the 
probabilities 

d-i 

P(Ai = Bj +n) = y~] P(Ai = m, B. } =n + m mod d) 

m=0 
d-1 

P(Ai ^ Bj+n) = V" P(Ai = m, Bj =l + m mod d) 

Ijtn 

(21) 

and to search for Bell inequalities written as linear com- 
binations of the P(Ai = Bj + n). This reduces consid- 
erably the number of Bell inequalities among which one 
must search in order to find the optimal one. The sec- 
ond idea is to search for the logical contradictions which 
force the Bell inequality to take a small value in the case 
of lhv theories. Thus the Bell inequality will contain 
terms with different weights, positive and negative, but 
the lhv theory cannot satisfy all the relations with the 
large positive weights. Once we had identified a candi- 
date Bell inequality, we ran a computer program that 
enumerated all the deterministic classical strategies and 
computed the maximum value of the Bell inequality. The 
deterministic classical strategies are those for which the 
probabilities PK 1 ...K Na L 1 ...L Nb are equal either to or to 
1 . In order to find the maximum classical value of a Bell 
expression, it suffices to consider them since the other 
strategies are obtained as convex combinations of the de- 
terministic ones. 

However when the number of settings, N a and Nb, 
and the dimensionality d increase, it becomes more and 
more difficult to find the optimal Bell inequalities using 
the above analytical approach. We therefore developped 
an alternative method based on the numerical algorithm 
which is used to find the threshold detection efficiency. 

The idea of this numerical approach is based on the 
fact that the probabilities for which there exists a solu- 
tion pkl to eqs. (|^,||||) form a convex polytope whose 
vertices are the deterministic strategies. The facets of 
this polytope are hyperplanes of dimension D — 1 where 
D is the dimension of the space in which lies the poly- 
tope (D is lower than the dimension (d + \) N ^+ N '> of 
the total space of probabilities due to constraints such as 
the normalizations conditions (Q) and the no-signalling 
conditions (Q)). These hyperplanes of dimension D — 1 
correspond to Bell inequalities. 

At the threshold 77*, the quantum probability 
belongs to the boundary, i.e to one of the faces, of the 



polytope determined by eqs The solution p^- L 

to these equations at the threshold is computed by our 
algorithm and it corresponds to the convex combinations 
of deterministic strategies that reproduce the quantum 
correlations. From this solution it is then possible to 
construct a Bell inequality. Indeed, the face F to which 
Px„ m belongs is the plane passing through the determin- 
istic strategies involved in the convex combination p Jj- L . 
Either, this face F is a facet, i.e. an hyperplane of di- 
mension D — 1, or F is of dimension lower than D — 1. 
In the first case, the hyperplane F correspond to the Bell 
inequality we are looking. In the second case, there is an 
infinity of hyperplanes of dimension D — 1 passing by F, 
indeed every vector v belonging to the space orthogonal 
to the face F determines such an hyperplane. To se- 
lect one of these hyperplanes liyng outside the polytope, 
and thus corresponding effectively to a Bell inequality, 
we took as vector v the component normal to F of the 
vector which connects the center of the polytope and the 
quantum probabilities when rj = 1 : Px„—i ■ Though this 
choice of v is arbitrary, it yields Bell inequalities which 
preserve the symmetry of the probabilities pQ M . 

As in the analytical method given above, we have ver- 
ified by enumeration of the deterministic strategies that 
this hyperplane is indeed a Bell inequality (ie. that it 
lies on one side of the polytope) and that it yields the 
threshold detection efficiency 77*. 



III. RESULTS 

Our results are summarized in table I. We now describe 
them in more detail. 



A. Arbitrary dimension, two settings on each side 

(N a =N b = 2). 

For dimensions up to 7, we found numerically that 
77* =1 = Tfl . The optimal measurements we found are 
identical to those maximizing the generalization of the 
CHSH inequality to higher dimensional systems Q , thus 
confirming their optimality not only for the resistance to 
noise but also for the resistance to inefficient detectors. 
Our values of 77* are identical to those given in [0 where 
r]* =1 has been calculated for these particular settings for 
2 < d < 16. 

We now derive a Bell inequality that reproduces an- 
alytically these numerical results (which has also been 
derived by N. Gisin |D|). Our Bell inequality is based 
on the generalization of the CHSH inequality obained in 
H . We recall the form of the Bell expression used in this 
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inequality: 

[d/2]-l 

rd,2x2 _ ST^ 

TV / . 

fc=0 



1 - 



2k 
d-1 



[P{A X = B ± + k) + P(B 1 = A 2 + k + 1) 

+P{A 2 = B 2 + k) + P{B 2 = Ai + *)] 
[P(Ai=Bi-A-l)+P(Bi = A 2 -A!) 
+P(A 2 = B 2 -k-l) + P(B 2 = Ai — k — 1)] 

(22) 

For local theories, I!// 2x2 < 2 as shown in [|| where the 
value of P// 2x2 (QM) given by the optimal quantum mea- 
surements is also described. In order to take into account 
"no-result" outcomes we introduce the following inequal- 
ities: 



j<2,2x2 _ jd,2x2 



2 4-^ 



P(A Z = 0, Bj = 0) < 2 (23) 



Let us prove that the maximal allowed value of / d > 2x2 for 
local theories is 2. To this end it suffices to enumerate 
all the deterministic strategies. First, if all the local vari- 
ables correspond to a "result" outcome then I!// 2x2 < 2 



and /, 



d,2x2 



B; 



so that 



jd,2x2 <- 2; if on e of the local variables is equal to 
then again I^ 2x2 < 2 (since the maximal weight of a 
probability in l^.: 2x2 is one and they are only two such 



d,2x2 
00 



probabilities different from zero) and I, 
are two outcomes, then I^- 2x2 < 1 and I, 
while if there are three or four 

t<2, 2 x 2 ,~ (-\ 



0; if there 



(2.2x2 



< 1; 



then lf- 2x2 = and 



eq. (| 



Note that the inequality ( |23| ) obeys the condition 
= c, hence it will provide a bound on r]^ x . Using 
, we obtain the value of rj% x : 



I?f x2 (QM) 



(24) 



Inserting the optimal values of I^ 2x2 (QM) given in Jq] 

this reproduces our numerical results and those of & 

As an example, for dimension 3, I 3 ]? x2 (QM) = 2.873 so 

that r% x = 0.8209. When d -> oo, (ph gives the limit 
v % - - - . - 1 — * 



- VA - 0.8049. 



B. 



3 dimensions, 2x3 settings. 



We have derived a Bell expression associated to these 
measurements: 



t3 : 2x3 



\P{M - 

+P{A 2 = B l 
- [P{A\ ^ Bi) - 
+P(A 2 ^ B l 



B x ) + P(A 1 = B 2 ) + P(A 1 = B 3 ) 
- 1) + P(A 2 = B 2 + 2) + P(A 2 = B 3 )} 

P{A 1 ± B 2 ) + P(A 1 ± B 3 ) 
- 1) + P(A 2 ^B 2 + 2) + P(A 2 ? B 3 )] 

(25) 



The maximal value of I^ 2x3 for classical theories is 2 
since for any choice of local variables 4 relations with a + 
can be satisfied but then two with a - are also satisfied. 
For example we can satisfy the first four relations but 
this implies A 2 = B 2 + 1 and A 2 = B 3 + 1 which gives 2 
minus terms. The maximal value of I^ 2x3 for quantum 
mechanics is given for the settings described above and 
is equal to I 3 f x3 (QM) = 10/3. To take into account 
detection inefficiencies consider the following inequality: 



r 3,2x3 



r 3,2x3 



r 3,2x3 



/i 2x3 < 2 



where 



and 



r 3,2x3 



r 3,2x3 



,B< 



(26) 



(27) 



(28) 



(J r is taken equal to zero). The principle used to show 
that / 3 < 2x3 < 2, is the same as the one used to prove 
that I d - 2x2 < 2. For example if A 1 = then I 3 < 2x3 < 3, 



r3 : 2x3 



= — 1 and I, 



3,2x3 



= so that 7 3 < 2x3 < 3 - 1 = 2. 

vi i 1 vv i — -n 

From ( |26[ ) and the joint probabilities ( |20| ) for the optimal 
quantum measurements we deduce: 



G 



■'/* 
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3 



— ~ 0.8182 



(29) 



in agreement with our numerical result. 

Note that in |T(| , an inequality formally idendical to 
( p5j ) has been introduced. However, the measurement 
scenario involve two measurements on Alice's side and 
nine binary measurements on Bob's side. By grouping 
appropiately the outcomes, this measurements scenario 
can be associated to an inequality formally identical to 
(25) for which the violation reaches 2^3. According to 
( 29 ) , this result in a detection efficiency threshold rj^ x of 
6/(2^3 + 4) w 0.8038. 



For three-dimensional systems, we found that adding 
one setting to one of the party decreases both -q x=1 and 
rji x from 0.8209 to 0.8182 (In the case of d = 2, it is 
necessary to take three settings on each side to get an 
improvement). The optimal settings involved are 4>a 1 = 
(0,0,0), cb A2 = (0,2^/3,0), <j> Bx = (O.tt/3,0), 0b 2 = 
(0, 2tt/3, -tt/3), <j>B 3 = (0, -it/3, -it/3). 



C. 3 settings for both parties 

For 3 settings per party, things become more surpris- 
ing. We have found measurements that lower rj x=1 and 
rj^ x with respect to 2 x 2 or 2 x 3 settings. But con- 
trary to the previous situations, rj x=1 is not equal to r/^ x , 
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and the two optimal values are obtained for two different 
sets of measurements. We present in this section the two 
Bell inequalities associated to each of these situations for 
the qubit case. Let us first begin with the inequality for 



J£ 3x3 - A = E(A 1 ,B 2 ) + E(A 1 ,B 3 ) + E{A 2 ,B 1 ) 
+E(A 3 , Bi) - E{A 2 , B 3 ) - E{A 3l B 2 ) 

-\p{Ai ± B x ) - ^P(A 2 + B 2 ) - ^P(A 3 ^B 3 )<2 

(30) 

where E(A U B ) = P(A l = B 3 ) - P(A l ^ Bj). As 
usually, the fact that J r 2 ; 3x3 < 2 follows from consid- 
ering all deterministic classical strategies. The maxi- 
mal quantum mechanical violation for this inequality is 
3 and is obtained by performing the same measurements 
on both sides A\ = Bi, A 2 — B 2 , A 3 — B 3 defined 
by the following phases: (j>A 1 — (0,0), 4>a 2 = (0, tt/3), 
4>a 3 = (0, — vr/3). It is interesting to note that this in- 
equality and these settings are related to those considered 
by Bell [Q and Wigner jl?]] in the first works on quantum 
non-locality. But whereas in these works it was necessary 
to suppose that Ai and Bj are perfectly (anti-)correlated 
when i = j in order to derive a contradiction with lhv 
theories, here imperfect correlations P(Ai ^ Bi) > can 
also lead to a contradiction since they are included in the 
Bell inequality. 

If we now consider "no-result" outcomes, we can use 
Irr 3x3 ' X without adding extra terms and the quantum 
correlations obtained from the optimal measurements vi- 
olate the inequality if 



The particular form of the terms /0^ 3x3 ' VA , 



r 2,3x3,VA , 
1 a and 

r0 



Inn* 3 ' A is given in the Appendix. The important point 



00 _ 
isthatJ\ Jife 



K+c. 



(0), (|) and (§g), we thus deduce 

4 



8/3 and £ 4J c 00 = 2. From 



\/A 



r 2,3x3,VA 



{QM) 



(34) 



The measurements that optimize the former inequality 
( |30| ) give the threshold r)^ x = 16/19. However these 
measurements are not the optimal ones for (|32|). The 
optimal phase settings are given in the Appendix. Using 
these settings it follows that I^f x3 ' yx (QM) = 3.157 and 
»# A ~ 0.8217. 

One may argue that the situation we have presented 
here is artificial and results from the fact that we failed 
to find the optimal inequality valid for all lambda which 
would otherwise have given 
identical to the threshold rj x ' 

the case since for A > 1 and rj > r) x=1 there exists a lhv 
model that reproduces the quantum correlations. This 
lhv model is simply given by the result of the first part 
of our algorithm described in II D. 



a threshold ^ A = 0.8165 
1 . However, this cannot be 

A=l 



D. More settings and more dimensions 

Our numerical algorithm has also yielded further im- 
provements when the number of settings increases or the 
dimension increases. These results are summarized in 
Table I. For more details, see the Appendix. 



W > 



r 2,3x3,A 



(QM) 



(31) 



IV. CONCLUSION 



Taking A = 1, we obtain r] X=1 = ^flfi ~ 0.8165. For 
smaller value of A, r\ x increase until r\ x = 16/19 is reached 
for A ~ 0.9401. At that point the contradiction with local 
theories ceases to depend on the production rate A. It is 
then advantageous to use the following inequality 

7 2,3x3,va = l E (A 1 ,B 2 ) + ^E(A U B 3 ) + ^E(A 2 ,B 1 ) 

+ ^E(A 3 ,B 1 ) - ^E(A 2 ,B 3 ) - ^E(A 3 ,B 2 ) 

* Si) - ^P(A 2 + B 2 ) - ^P(A 3 ^B 3 )<2 

(32) 

This inequality is similar to the former one (^0|) but the 
symmetry between the E(Ai, Bj) terms has been broken: 
half of the terms have an additional weight of 1/3 and 
the others of —1/3. The total inequality involving "no- 
result" outcomes is 

r2.3x3.VA _ t-2, 3x3, VA , r2.3x3.VA , r2.3x3.VA , r2.3x3.VA < 9 

(33) 



In summary we have obtained using both numerical 
and analytical techniques a large number of Bell inequal- 
ities and optimal quantum measurements that exhibit an 
enhanced resistance to detector inefficiency. This should 
be contrasted with the work (reported in [||, |j| ) devoted 
to searching for Bell inequalities and measurements with 
increased resistance to noise. In this case only a single 
family has been found involving two settings on each side 
despite extensive numerical searches. Thus the structure 
of Bell inequalities resistant to inefficient detectors seems 
much richer. It would be interesting to understand the 
reason for such additional structure and clarify the origin 
of these inequalities. 

It should be noted that for the Bell inequalities we 
have found, the amount by which the theshold detector 
efficiency 77* decreases is very small, of the order of 4%. 
This is tantalizing because we know that for sufficiently 
large dimension and sufficiently large number of settings, 
the detector efficiency threshold decreases exponentially. 
To increase further the resistance to inefficent detector, 
it would perhaps be necessary to consider more general 
measurements than the one we considered in this work or 



9 



use non-maximally entangled states (for instance, Eber- 
hard has shown that for two-dimensional systems, the 
efficiency threshold 77* can be lowered to 2/3 using non- 
maximally entangled states fL8[ ). There may thus be a 
Bell inequality of real practical importance for closing the 
detection loophole just behind the corner. 
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d = 2, N A = 3, N B = 3, A 

Bell inequality: 

I 2 < 3 * 3 > x = E(Ax,B 2 ) + E(A lt B 3 ) 

+ E(A 2 ,B 1 ) + E(A 3 , Bx) - E{A 2l B 3 ) 

-E(A 3 ,B 2 )-^P(A 1 ^B 1 ) 

- ^P(A 2 B 2 ) - ^P(A 3 ^B 3 )<2 

where E{A i ,B j ) = P(A, = Bj) - P(A, ^ Bj). 
Optimal phase settings: 

<t> Al = (0,0) 4>a 2 = (O.tt/3) 4>a 3 = (0,-tt/3) 

<t> Bl = (0,0) <j) B2 = (0,7T/3) <t>B 3 = (0.-7T/3) 



APPENDIX A 

For completeness, we present here in details all the Bell 
inequalities and optimal phase settings we have found. 
This includes also the results of Table I which have not 
been discussed in the text. 

• N A = 2, N B = 2, VA 

Bell inequality: 



[d/2]-l 



T d,2x2 



E 1 



fc=0 



2k 
d-l 



\P{Ai = Bi + k) + P(Bi = A 2 + k + 1) 

+ P(A 2 = B 2 + k) + P{B 2 =Ai + k)] 
- [P(Ai =Bi-k-l) + P(Bi =A 2 -k) 

+ P(A 2 = B 2 -k-l) + P(B 2 = A x - k - 1)] 

1 2 

+ 5 E P(A i = 0,fl i =0)<2 



Maximal violation: I 2 < 3X6 > A (QM) = 3 
Detection threshold: — \/ T\ 



• d = 2, N A = 3, N B = 3, VA 

Bell inequality: 



r 2.3x3.VA 



= -E(A 1 ,B 2 ) + -E(A 1 ,B 3 ) 



3 " " 3 

+ ^E(A 2 , Bx) + ^E(A 3 ,Bx) - ^E(A 2 ,B 3 ) 

-^E(A 3 ,B 2 )-^P(Ax^Bx) 

- ^P(A 2 ? B 2 ) - ^P(A 3 ? B 3 ) 

- ^F % (Ax,B 2 ) - ^F % (A 2 ,B 3 ) - ^F (A 3 ,Bx) 
+ ^Fu(A 3 ,B 2 ) + ^P(A 2 =%,Bx^ 0) 

+ |p(A x + 0, B 3 = 0) + ^P(Ax =®,Bx = 0) 
+ \p{M = 0, Bx = 0) + |i>(Ai = 0, B 3 = 0) < 2 



Optimal phase settings: 



<Mi) = o <Mi) = fj 
toiV) = jaJ 4>B 2 {i) = -jaJ 



Maximal violation: 



[d/2]-l 



I d ' 2x2 (QM) =4d ^ ( 1 -jh) {qk - q - k - l) 



k=0 



2k 



where E{A l ,B j ) = P(A t = Bj) - P{A % ^ Bj) and 
F$(Ai,Bj) = P{A t = Q,Bj + 0) + P(A t ? <D,Bj = 
0)+P(A i = 0,B j =0). 

Optimal phase settings: 

4> Al = (0,0) tf, Aa = (0, 1.3934) 
<Pa 3 = (0, -0.7558) 

4> Bl = (0, 0.5525) (j)B 2 = (0, 1.3083) 
4> Ba = (0, -0.8410) 



where q k = 1/ (2d 3 sin 2 [7r(A: + l/4)/d]). 
Detection threshold: 77^ A = Jd ,2x2(q A/ ) +2 



Maximal violation: 7 2 < 3x3 ' VA (QM) = 3.157 
Detection threshold: r^ x = 0.8217 
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• d = 2, N A = 3, N B = 4, VA 

Bell inequality: 

/2 ,3X4,VA = _ p{Ai ^ _ p{Ai _ p{Ai ^ 

+ P(A 2 = Bi) + P(A 2 = B 2 ) - P(A 2 ? B 3 ) 
+ P(A 2 ? B^ - P(A 3 = B x ) + P{A 3 = B 2 ) 

- P(A 3 ? B 2 ) + P(A 3 ? B 3 ) - P(A 3 = Bi) 
+ P{A l £ 0, Si = 0) + P(A 2 = 0, B l ± 0) 

- P(A 3 £ 0, Bt = 0) - P{Ax = 0, B 2 ^ 0) 

+ P(Ai = 0, B x = 0) + P{A 2 = 0, B 2 = 0) < 2 



Optimal phase settings: 

^ = (0,0) 0a 2 = (0,0.7388) 
<Pa 3 = (0,2.1334) 

4> Bl = (0,-0.1347) 4>b 2 = (0,1.2938) 

<j) B3 = (0, -0.0757) B4 = (0, -1.0891) 



Maximal violation: I 2 ' 3xi (QM) = 2.8683 
Detection threshold: r?* VA = 0.8216 

• d = 2, N A = 4, N B = 4, VA 

Bell inequality: 

/2 ,4X4,VA = _ p(Ai = + p{Ai _ p{A2 = Bi) 

- P(A 2 = B 2 ) + P(A 2 ? Bi) + P(A 3 ± Si) 

- P(A 3 ^ B 2 ) - P(A 3 + B 3 ) - P(Ai ± Bi) 

- P(Ai = B 2 ) - P(Ai = B 3 ) + P(Ai + Bi) 
+ P{A l £ 0, Bi = 0) - P(A 4 ^ 0, B l = 0) 

+ P{A X = 0, Bi = 0) + P(Ai = 0, Bi = 0) < 2 



Optimal phase settings: 

^=(0,0) <j>A 2 = (0,0.0958) 

<f>A 3 = (0,2.1856) <f> Ai = (0,4.5944) 

4> Bl = (0,4.0339) <j>B 2 = (0,3.3011) 

0s 3 = (0, 2.2493) B4 = (0, 2.3454) 



Maximal violation: / 2 < 4x4 (QM) = 2.8697 
Detection threshold: r^ A = 0.8214 



• d = 3, N A = 2, N B = 3, VA 

Bell inequality: 

/3 ,2x 3 ,a = + [p(Aj = Bi) + P(Ai = B 2 ) + P(A X = B 3 ) 
+ P(A 2 = B 1 + 1) + P{A 2 = B 2 + 2) + P(A 2 = B 3 )} 
- [P(Ai ? Bi) + P(A 1 ± B 2 ) + P(A 1 + B 3 ) 
+ P(A 2 £ B l + 1) + P(A 2 ? B 2 + 2) + P(A 2 ? B 3 )] 

-i]TP(^=0,B,^0) 
+ i^P( J 4 l =0,B,=0)<2 



Optimal phase settings: 

cj> Al =(0,0,0) </u 2 = (0,2tt/3,0) 

Bl = (0, tt/3, 0) 0b 2 = (0, 2tt/3, -tt/3) 

<^>b 3 = (0, -tt/3, -tt/3) 



Maximal violation: 7 3 < 2x3 (QM) = f 
Detection threshold: r^ A = £ ~ 0.8182 



• d = 3, N A = 3, AT B = 3, A 

Bell inequality: 

/ 3 ' 3x3 ' A =Bi(Ai,B 2 ) + B 2 (^i,B 3 ) 
+ B 2 (A 2 , Bi) - B 2 (A2, B 3 ) + E l (A 3 , Bi ) 
-E 1 {A 3 ,B 2 )-P{A 1 ^B 1 ) 
- P(A 2 £ B 2 ) - P(A 3 ? B 3 ) < 2 



Optimal phase settings: 

<P Al = (0, 0, 0) <f> A2 = (0, 2tt/9, 4tt/9) 

cj> As = (0, -27T/9, -47T/9) 

Bl = (0, 0, 0) 0s 2 = (0, 2tt/9, 4tt/9) 
<Pb 3 = (0, -2tt/9, -4tt/9) 



Maximal violation: I 3 < 3x3 (QM) = 3.0642 
Detection threshold: 7/ A = 3 g 42A 
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• d = 3, N A = 3, N B = 3, VA 

Bell inequality: 

/3 ,3x 3 ,va = -l P ( Al = Pi) - |p(Ai = Pi + 2) 
+ P(A 1 = B 2 ) + |p(Ai = B 2 + 1) - |P(^ = B 3 ) 

- P(Ai = P 3 + 2) + ^P(A 2 = Pi) - 2P(A 2 = Pi + 1) 

- ^P(A 2 = P 2 ) + 2P(.4 2 = P 2 + 1) - P(A 2 = B 3 + 1) 

- ^P(A 2 = P 3 + 2) - yP(^3 = Pi) - 2P(A 3 = Pi + 2) 
+ ^P(A 3 = B 2 ) + 2P{A 3 = P 2 + 1) + ^P(A 3 = P 3 ) 
+ P(A 3 = P 3 + 2) + |p(^! ^ 0, Pi = 0) 

- jjp(A 2 ^ 0, Pi = 0) - 2P(A 3 ^ 0, Pi = 0) 

+ 2P(A 1 £ 0, P 2 = 0) + |p(A a = 0, Pi = 0) 
+ 2P{A 1 = 0,P 2 = 0) < 11/3 

Optimal phase settings: 

4> Al = (0, 0, 0) (j>A 2 = (0, 1.4376, 2.8753) 
4> A3 = (0,0.5063,1.0125) 

4> Bl = (0,2.0452,4.0904) <j> B2 = (0,2.9758,-0.3315) 
4> Ba = (0,1.3839,2.7678) 

Maximal violation: I 3 ' 3x3 (QM) = 5.3358 
Detection threshold: rj* x = 0.8146 

• d = 4, N A = 2, N B = 3, VA 

Bell inequality: 
/4 ,2x3,va = p^ Ai =Bl + i) + 2P(A X = Pi + 2) 
+ 2P{A X = P 2 ) + P{ Al = P 2 + 1) + 2P{A l = P 3 ) 

+ 2P(^ 2 = Pi + 1) + P(A 2 = Pi + 2) + P(A 2 = P 2 ) 

+ 2P(A 2 = P 2 + 1) + 2P(^1 2 = P 3 + 2) 
+ ^P(A 4 = 0,Pi^0) + i^P(A t = 0,P 2 ^0) 

i i 

+ i^P(A 4 = 0,P 3 ^0) + ^P(Ai^0,P 1 = 0) 

i i 

+ i^P(A 2 ^0,Pi =0) + ^P(Ai = 0,Pi = 0) 

i 

+ |p(Ai = 0, P 2 = 0) + |p(^ = 0, P 3 = 0) 
+ ^P(A 2 = 0, Pi = 0) + ip(A 2 = 0, P 2 = 0) 
+ ip(^ 2 = 0,P 3 = 0)<8 



Optimal phase settings: 
<f> Al = (0,0,0,0) 

<^4 2 = (0,-1.1397,2.0019,3.1416) 
<f> Bl = (0,1.7863,-0.5698,2.3562) 
ct) B2 = (0,0.2155,5.7133,0.7854) 
(f> B3 = (0,1.0009,1.0009,0) 



Maximal violation: I 4 ' 2x3 (QM) = 9.4142 
Detection threshold: r% x = 0.8093 



• d = 4, N A = 3, N B = 3, VA 

Bell inequality: 

7 4,3X3,VA = _ P ( Al =B 1 + 2) + P{A X = Pi + 3) 

+ 2P{A 1 = B 2 + 1) - P(Ai = P 2 + 2) - P(Ai = P 3 ) 

- 3P(Ai = P 3 + 1) - 2P(Ai = P 3 + 2) - P(A 2 = Pi) 

+ P(A 2 = Pi + 1) - P(A 2 = P 2 + 1) + P(A 2 = P 2 + 2) 

+ 2P(A 2 = P 3 + 3) + 2P(A 3 = Pi + 1) + P(A 3 = B 2 ) 

- 2P(A 3 = P 2 + 2) - P(A 3 = P 2 + 3) + 2P(A 3 - P 3 ) 

+ P(A 3 = P 3 + 2) + ^ P(Ai = 0, Pi ^ 0) 

+ P(A! ? 0, Pi = 0) + P(Ai ^ 0, P 2 = 0) 

- P{A X = 0, P 3 ^ 0) + P(A 3 = 0, Ps ^ 0) 
+ P(A 3 £ 0, P 3 = 0) + 2P(Ai = 0, Pi = 0) 
+ P{A X = 0, P 2 = 0) + P(A 2 = 0, Pi = 0) 

+ P(A 3 = 0, Pi = 0) + P(A 3 = 0, P 3 = 0) < 6 



Optimal phase settings: 
<p Al =(0,0,0,0) 

0a 2 = (0,-1.2238,-1.1546,3.9048) 
0a 3 = (0,3.1572,3.8330,0.7070) 
(j) Bi = (0,-0.9042,1.7066,0.8025) 
(j) B2 = (0,2.5844,3.6937,-0.0051) 
<j) Ba = (0,4.1396,3.0022,7.1419) 



Maximal violation: I 4 > 3x3 (QM) = 7.5576 
Detection threshold: = 0.7939 
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